Characterizing Potentials by a Generalized Boltzmann Factor 
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Based on the concept of a nonequilibrium steady state, we present a novel method to experi- 
mentally determine energy landscapes acting on colloidal systems. By measuring the stationary 
probability distribution and the current in the system, we explore potential landscapes with barri- 
ers up to several hundred fca T. As an illustration, we use this approach to measure the effective 
diffusion coefficient of a colloidal particle moving in a tilted potential. 
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Introduction. - The interaction of soft matter systems 
with potential landscapes created by optical tweezers 
plays a key role for, e.g., mechanical flexibility measure- 
ments of single biomolecules or molecular motors [l|, Q , 
guiding of neuronal cells [|| , or phase transitions of col- 
loidal monolayers on patterned substrates [1, Q . In ad- 
dition, extended optical lattices can be used as sorters 
for microscopic particles [|| or as microoptomechanical 
devices such as Couette rheometers Q. Currently, no 
theories are available which can be used to directly cal- 
culate optical trapping forces on macromolecules. Thus 
the precise calibration of optical forces is a central issue 
in many experiments. 

The simplest method to determine an optical potential 
V(r) is to measure the equilibrium distribution p eq (r) of 
a highly diluted colloidal system at position r. From the 
inverted Boltzmann factor 



V(r) = -feTlnp eq (r) 



(1) 



one directly obtains the underlying potential V(r) with 
a typical energy resolution on the order of 0.1 T [3, Hi, 
where fce is Boltzmann's constant and T the temperature 
of the surrounding fluid. This technique, however, is only 
applicable to potential depths up to ~ 7 ks T which are 
effectively sampled by Brownian particles in equilibrium. 
For larger trapping potentials, optical forces are typically 
calibrated indirectly by taking advantage of Stokes law 
which relates the particle velocity to the friction force 
exerted by the surrounding solvent molecules. Accord- 
ingly, from the drift velocity of a particle, the under- 
lying potential can be reconstructed SO, EH. Al- 
ternatively, within the drag force method, V{r) can be 
determined from the particle's displacement upon mov- 
ing the sample stage (and thus the liquid) with known 
velocity P, 0, [H, El- In contrast to Eq. (JTJ, however, 
the latter two nonequilibrium methods neglect thermal 
fluctuations since only mean values of particle velocities 
or displacements are considered. While such fluctuations 
can be neglected at large trapping forces, this is no longer 
justified for external forces with strengths comparable to 
those exerted by fluctuating Brownian forces. 

In this paper, we introduce a potential reconstruction 
method based on a generalization of Eq. ([T]) to nonequi- 
librium conditions. This is experimentally realized by 



generating a nonequilibrium steady state (NESS) for a 
colloidal particle in a one-dimensional (toroidal) poten- 
tial landscape. By measuring the stationary probability 
distribution and the current in the system, we can reli- 
ably calibrate potentials wells between a few tens up to 
several hundreds of ks T. 

Potential reconstruction. - Our method is based on a 
generalization of the Boltzmann factor inversion (JTJ) to 
nonequilibrium. The effectively one-dimensional motion 
of the particle along a toroidal trap is governed by a 
Langevin equation 



r 1 F(x)+at) 



(2) 



with x the spatial coordinate and Q{t) representing ther- 
mal noise with correlations (C(i)C(O) = 2(feT H)&{t — 
t'), where 7 is the friction coefficient. The force F(x) = 
—V'(x) + f exerted on the particle stems from two 
sources, the gradient of the periodic potential V(x + L) = 
V(x) and a nonconservative driving force /. 

We define a pseudo-" potential" 4>(x) by writing the 
nonequilibrium steady state probability distribution as 
Ps(x) = cxp[— 4>(x)] resembling the Boltzmann factor. 
The stationary probability current through the toroid is 
given as 



j s = 1 - 1 {F(x) Ps (x)-k B T P ' s (x)} 



(3) 



which is constant in one dimension. We introduce the 
local mean velocity 



v s (x) = j s /p s (x) 



and obtain uA 



1 v s {x) = F{x)+k B T(j>'{x) = -V'(x)+f+k B T^{x). (4) 
Integration of Eq. @ leads to the potential 

V{x) = k B T<j>{x)+ [ [f- 1Vs (x)]dx (5) 
Jo 

up to an irrelevant additive constant. Using the defini- 
tions of 4>(x) and v s (x) we finally arrive at 



V{x) = -kBTlnp s (x) + fx - jj s / p s (x)dx. (6) 
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FIG. 1: (color online) Stationary probability distribution 
p s (x) (black bars) and the pseudopotential 4>{x) (dashed line), 
measured at a driving force of 34 ka T/fim, pointing along the 
negative x-direction. The potential V(x) (dotted line) is de- 
termined according to Eq. ([6|. Note that 4>(x)is multiplied by 
a factor of 10 and shifted vertically to enhance visualization. 

Hence, the stationary probability p s (x) and the local 
mean velocity v s (x) determine the potential V(x). The 
driving force / can be determined by setting x = L in 
Eq. j5]) and using the periodicity of the potentials V(x) 
and 4>{x) as 

f=lJ L v s {x)dx. (7) 

In thermal equilibrium both j s and / vanish and Eq. 
reduces to the inverted Boltzmann factor Eq. ([1]) . There- 
fore Eq. (JHJ) can be understood as an extension of the 
Boltzmann factor to noncquilibrium stationary states. 

Experiment. - For an experiment exploiting Eq. ([5]), 
we use a scanning optical tweezers setup as described 
in detail elsewhere [15(- A laser beam (A ~ 532 nm) 
is deflected on a pair of galvanomctric mirrors and fo- 
cused with a lOOx, NA=1.3 oil immersion objective from 
below onto a silica bead immersed in water (diameter 
d ~ 1.85 fim). Upon periodic modulation of the angular 
mirror positions we obtain a three dimensional toroidal 
laser trap with a torus radius of R ~ 3.95 fim. At our 
driving frequencies ~ 100 Hz, the particle can not fol- 
low directly the rotating laser trap. Instead, every time 
the particle is passed by the laser tweezers, it experiences 
a minute kick along the rotation direction whose strength 
depends on the laser intensity Iq [HI- Because the par- 
ticles trajectory is monitored with video microscopy at 
a sampling rate of 20 Hz, single kicking events are not 
resolved and the driving force / alon g th e angular di- 
rection can be considered as constant [15j . In addition, 
the intensity of the laser is weakly modulated along the 
toroidal trap. This is achieved with an electro-optical 
modulator (EOM) being controlled by a function gener- 



FIG. 2: (color online) Averaged potential (solid points). The 
grey curve is a fit of Eq. © to the averaged potential. Inset: 
Averaged potentials obtained for driving forces / (= 34, 43, 
57 and 73fc B T//im). 

ator which is synchronized with the scanning motion of 
the mirrors. This intensity modulation I m (x) leads to 
a periodic potential V(x) with x the arc-length coordi- 
nate along the circumference of the torus. It has been 
demonstrated that the resulting optical forces of such an 
intensity modulated scanned laser tweezers exerted on a 
colloidal particle correspond to those of a tilted periodic 
potential [l5[ . 

To experimentally demonstrate that V(x) can be ob- 
tained under non-equilibrium steady state conditions, the 
intensity of the scanned laser tweezer along the toroidal 
trap was varied according to 

I(x) = I +I m am(n— ). (n = 4) (8) 

Fig. [1] shows the steady state probability distribution 
p a {x) as obtained from the particles trajectory and the 
corresponding pseudo-potential 4>(x) for Iq = 44 mW 
and I m = 10 mW. Together with the driving force 
/ = 34fceT//Ltm as determined from the measured lo- 
cal mean velocity v s (x) (cf. Eq. {?])) we finally arrive 
at the potential V(x) which is also plotted in Fig. [T] as 
dotted line. Clearly, under NESS conditions the minima 
and maxima of <fi(x) and V(x) do not coincide. In ad- 
dition, <fi(x) varies in a less pronounced way than V(x) 
because p s {x) is broader than it would be in equilibrium. 
On top of the intensity modulation according to Eq. ([5]) 
we observe a constant, small variation of the potential 
with 27r-periodicity caused by minute optical distortions 
in our setup. Since we are only interested in the local 
shape of the potential, in the following we only consider 
potentials, where V(x) is averaged over the four exter- 
nally applied periods. The averaged potential, which is 
plotted as solid bullets in Fig. [2l is in excellent agreement 
with a sinusoidal fit to 

V(x) = ^ S m(A (9) 
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FIG. 3: (color online) Potential depth as a function of the 
modulation amplitude I m . The different symbols correspond 
to data acquired at different driving forces. (B:43, A:57 and 
T:73feT//xm) 



as theoretically expected for the optical potential in case 
of sinusoidal intensity variations [161 ] . 

To demonstrate the robustness of our approach in char- 
acterizing equilibrium potentials under NESS conditions, 
we systematically varied the driving force f — 34, 43, 57 
and 73 ksT / fj.ni while keeping V(x) unchanged. Exper- 
imentally, this is achieved by changing Iq with all other 
parameters in Eq. ([8]) fixed. The measured potentials 
plotted in the inset of Fig. [5] clearly fall on top of each 
other and thus demonstrate that the measured V(x) is 
independent of /. 

Similarly, the potential amplitude Vq can be changed 
by variation of I m . Fig. [3] shows the potential depth as 
a function of I m , measured at different driving forces / 
(marked by different symbols). As expected, we find a 
linear dependence between Vq and I m independent of /. 
This shows that reconstruction of equilibrium potentials 
can be reliably performed for a wide range of driving 
forces. For practical purposes, however, the driving force 
should not exceed a certain range where the lower limit 
is reached when the particle only rarely surmounts adja- 
cent potential barriers and thus cannot sample the entire 
landscape. For very large /, the probability distribution 
becomes rather flat and very long sampling times are re- 
quired to accurately measure V(x). 

Diffusion in tilted periodic potentials. - Having demon- 
strated the validity of our approach to reconstruct equi- 
librium potentials under NESS conditions, in the follow- 
ing we will exemplarily apply this method to the problem 
of giant diffusion. It has been shown theoretically [lTj 
and experimentally [HI, [HI that the effective diffusion co- 
efficient D c g of a Brownian particle moving in a tilted pe- 
riodic potential U (x) = V(x) + fx exhibits a pronounced 
maximum as a function of the driving force /. Until now 
experiments were not able to match quantitatively the 
theoretical predictions. With the ability to character- 



ize the underlying potential landscape in detail, we can 
quantitatively test the theoretical behavior of D c g. 

The effective diffusion coefficient is easily obtained 
from the particle trajectory according to 



D eS = lim D(t), D{t) 



2t 



(10) 



This expression takes into account both the thermal dif- 
fusion and the drift motion evoked by the tilt of the po- 
tential. It is therefore applicable to both equilibrium and 
nonequilibrium conditions. Depending on the strength of 
the driving force, three regimes can be distinguished: (i) 
At small /, the particle is largely confined to the poten- 
tial V(x). Thus D ff < Dq with Do the diffusion coef- 
ficient of a free particle, (ii) Around a critical force / c , 
a considerable enhancement of the thermal diffusion oc- 
curs, i.e. D c s > Dq [l?]]. (iii) In the limit of very large 
/ the potential becomes irrelevant and -D e ff eventually 
approaches Dq. 

Our results are shown in Fig. 01 where we have chosen 
the same sinusoidal potential as above (see Eq. [9]) with 
typical amplitudes between 10 — 20 k^T [2(J. Since the 
infinite time limit required to calculate D c h cannot be 
realized in experiments, we first plotted the right hand 
side of Eq. (|10p as a function of time to determine when 
this expression saturates. The inset of Fig. [4] shows the 
result obtained for / ~ 6.6kT//im and Vq ~ 10.3 feT. 
After an initial peak, the curve converges to the corre- 
sponding long-time value. A closer inspection reveals two 
damped oscillations whose periods are easily explained: 
The short oscillation time T\ ~ 1.7s corresponds to the 
mean residence time of the particle within one minimum 
while the other oscillation with r 2 ~ 6.6 s ~ 4ti equals 
the mean revolution time of the particle along the torus. 
After about t > 15 s, both oscillations have essentially 
decayed to the long-time value corresponding to D e ff . 

Fig. [4] shows the normalized effective diffusion coeffi- 
cients for potential depths of 10.3 and 14.4 k^T . Both 
curves show a peak clearly indicating the enhancement 
of thermal diffusion in tilted periodic potentials. With 
increasing potential strength we observe a shift of the 
D g(f) curve towards larger forces. The values of / c and 
Dcs(fc) sensitively depend on the shape of the poten- 
tial and are theoretically predicted as f c = 2Vq/R and 
D cS (fc) = 0.0696ZViU 7r 3 )( 2 / 3 ) 0. 

A comparison with our data is shown in Tab. [U While 
the predicted critical force f c is in rather good agree- 
ment with the experimental data, the theoretical values 
systematically underestimate D c g(f c )/Do by a factor of 
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TABLE I: Comparison between experimentally determined 
and theoretically predicted position and height of the giant 
diffusion peak. 
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those distortions, the good agreement for the critical tilt 
can be explained. 

Concluding perspective. - So far, we have demon- 
strated a novel method to reconstruct equilibrium po- 
tentials on the basis of the stationary probability distri- 
bution. In particular in one-dimensional NESS condi- 
tions, this quantity is easily determined experimentally, 
because the stationary current j s is constant. When the 
method is extended to higher dimensions in the pres- 
ence of nonconservative force fields f(r), in addition to 
the steady state probability p s (r) the local mean veloc- 
ity v s (r) = j s (r)/p s (r) is required. Experimentally, this 
quantity is obtained by averaging the velocity of parti- 
cles passing r. Then the actual potential could be re- 
constructed through integration along open paths C(r) 
starting at an arbitrary but fixed initial point and ending 
in r, leading to 



FIG. 4: (color online) Normalized diffusion coefficient 
D e g/Do vs external driving force. The data was obtained 
for potential depths Vb of 10.3 (►) and 14.4 feT (■). Inset: 
The function D(i) in Eq. (TjJjJ versus time t (/ ~ 6.6 kT/pm, 
Vo ~ 10.3 ksT). After 15s the long time limit is reached. 
The remaining small oscillations define the error in determin- 
ing D eff . 

about 0.7. The origin of this discrepancy is due to the 
aforementioned slight distortions along the toroidal trap 
which leads to local variations in the potential depth and 
thus affects the effective diffusion coefficient. Since at the 
same time the local shape of V{x) is hardly affected by 



V{v) = -k B T]np s {v) 

+ f [f(r( S ))- 7 v s (r( S ))]-rfr( S ). (11) 

JC(r) 

In summary, we have demonstrated a flexible method 
to characterize potentials using the generalization of the 
inverted Boltzmann factor. In contrast to equilibrium 
measurements, this allows to characterize laser potentials 
up to depths of several hundredth or even thousandths of 
&b T. Based on the determination of the stationary state 
probability distribution p s , this technique is easily appli- 
cable to different situations, e.g. topographical potentials 
and does not require fast data acquisition techniques. 
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